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ABSTRACT 



Our study begins by revisiting the photometric Butcher-Oemler effect using data 
from 7 CNOC clusters at 0.23<z<0.43. We construct the foreground/background cor- 
rected colour-magnitude diagrams for these clusters. Our analysis shows that the 
CNOC clusters reproduce the trend of increasing blue galaxy fraction with redshift as 
seen by Butcher & Oemler. We use these data to investigate the history of star for- 
mation in clusters by connecting these clusters as snapshots at different redshifts. We 
address two key issues. Firstly, we ask whether the simple fading and passive evolution 
of the blue galaxies is consistent with the properties of galaxies in nearby clusters, such 
as the Coma cluster. We find that the evolution of star forming field galaxies towards 
redder colour (once the star formation ceases on entry into the cluster environment) 
can successfully reconstruct colours and magnitudes of galaxies in the local cluster. 
There is no requirement for widespread disruption of these galaxies. Since the blue 
galaxies fade as they age, the fainter galaxies on the present-day colour-magnitude 
relation tend to have more extended star formation histories than their bright coun- 
terparts. However, this effect is not sufficient to cause a sizable age variation for the 
galaxies along the colour-magnitude relation, implying that the slope is dominated by 
variations in metal abundance. 

Secondly, we address the nature of the Butcher-Oemler effect itself. We compare 
the distribution of galaxies in the colour-magnitude diagrams and hence infer the 
evolution of the rate at which galaxies have arrived in the cluster. Models in which star 
formation is abruptly truncated as galaxies are accreted by the cluster have difficulty 
in reproducing the observed colour distribution. In contrast, if star formation declines 
on a 1 Gyr timescale after accretion, the galaxy accretion history we infer is consistent 
from cluster to cluster and matches well the distribution expected in simple theoretical 
models. The Butcher-Oemler effect is thus driven both by the declining star formation 
rates of field galaxies and by a decline in the rate at which fresh galaxies arc accreted 
by the cluster. 

Our study naturally leads to a comparison of the global star formation histories of 
galaxies in clusters and the field. We show that the star formation rate per galaxy mass 
for galaxies in cluster cores is significantly smaller than that of the field environment 
below z<l due to the truncation of star formation. However, the factor by which star 
formation is suppressed is dependent on the cluster accretion history. High quality 
observations of clusters at higher redshifts are needed to better define this relation. 

Key words: galaxies: clusters - galaxies: formation - galaxies: evolution - galaxies: 
stellar content 



1 INTRODUCTION 

It is 20 years since Butcher & Oemler (1978; 1984, here- 
after B084) first reported their startling results on the 



blue galaxy fractions of distant clusters. This and subse- 
quent photometric work suggested that the cores of distant 
clusters at redshift z>0.2 contain large numbers of blue 
galaxies (Couch & Newell 1984; Lubin et al. 1996; Rakos 
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& Schombert 1995), in stark contrast to the homogeneous 
red populations of local galaxy clusters (eg., Bower, Lucey 
& Ellis 1992; Terlevich, Bower & Caldwell 2000). There 
is an apparent contradiction between the narrow colour- 
magnitude (CM) relation of evolved galaxies in local clus- 
ters (which indicates a long period of passive evolution) and 
the blue colours of intermediate redshift galaxies (which in- 
dicates considerable, relatively recent star formation activ- 
ity). Bower, Kodama & Terlevich (1998, BKT98), Smail et 
al. (1998) and Balogh et al. (2000) have begun the process 
of modelling the combined dynamical and photometric evo- 
lution of clusters in order to reconstruct the galaxies' for- 
mation histories. We continue these themes in this paper. 

B084's original discovery of strong cluster evolution in- 
volved only the photometric properties of galaxies. We will 
refer to this aspect as the 'photometric Butcher-Oemler ef- 
fect'. Dressier & Gunn (1992) and Couch & Sharpies (1987) 
undertook spectroscopy initially with the motivation of con- 
firming spectroscopically that the blue galaxies were cluster 
members. They discovered that many of the blue cluster 
galaxies, and some red ones, had spectra of a type that 
are rare in local galaxy samples. These galaxies (referred to 
as E+A, PSG, a+k/k+a) have anomalously strong Balmer 
lines indicating an excess of A-stars superimposed on an K- 
star type spectrum. The galaxies are produced when star 
formation is abruptly truncated (Couch & Sharpies 1987). 
The strongest lined examples require a burst of star for- 
mation just before the truncation (eg., Barger et al. 1996; 
Poggianti & Barbaro 1996; Poggianti et al. 1999; Couch et 
al. 1998). We refer to the increase in the incidence of such 
spectra with redshift as the 'spectroscopic Butcher-Oemler 
effect' 

With the advent of HST, it has become possible to 
study the morphologies of galaxies in the distant clusters 
(Smail et al. 1997; Dressier et al. 1994; 1997; Couch et 
al. 1994; 1998). They found that the morphological mix in 
these clusters also evolved with redshift, creating a third 
variant of the Butcher-Oemler effect, which we will refer to 
as the 'morphological BO effect'. In particular, the fraction 
of SO galaxies in the cluster cores decreases rapidly with 
redshift (Dressier et al. 1997). Several attempts have been 
made to connect photometric/spectroscopic features and the 
morphology. Dressier et al. (1994) and Couch et al. (1994) 
showed that most of the blue B084 galaxies are actually 
normal spirals with the rest being interacting galaxies or 
mergers. Poggianti et al. (1999) compared the spectroscopic 
and the morphological effects, suggesting that the time scale 
of the morphological transition is longer than the A-type star 
evolution since the post star burst galaxies (k+a/a+k) are 
observed as normal spirals rather than SO galaxies. 

In this paper, we aim to clearly separate these three 
different (but related) evolutionary trends. We will focus on 
the photometric effect, returning to the questions raised by 
Butcher & Oemler's original papers. The advantage of the 
photometric approach is that it allows us to study large, 
complete samples of galaxies rather than a smaller subset of 
the population that is complicated by target selection and 
the strength of spectral features and often limited to brighter 
galaxies. Recent studies have tended to concentrate on par- 
ticular spectroscopic features to learn details about the na- 
ture of the star- burst /truncation phenomena seen in some 
cluster members (Barger et al. 1996; Poggianti et al. 1999; 



Balogh et al. 1999). On the other hand, we have made lit- 
tle progress on the evolution of the cluster population as a 
whole, as galaxies accrete from field, truncate their star for- 
mation, and fade and become red towards the CM relation 
ridge-line. A conspicuous contradiction exits between Pog- 
gianti et al. (1999) and Balogh et al. (1999) over the way star 
formation is suppressed in the cluster environment. While 
Poggianti et al. (1999) require an abrupt end to star forma- 
tion activity to explain the incidence of a+k/k+a spectra, 
Balogh et al. (1999) favor a slow decline in star formation in 
order to explain the radial dependence of the [On] emission 
line strength. 

To remedy this situation, we will look at the CM dia- 
grams of the distant clusters in more detail, examining how 
the whole cluster population evolves with time. The essential 
idea is that we determine the statistical star formation his- 
tories of galaxies in clusters by connecting a series of cluster 
snapshots at different redshifts. The process has been exten- 
sively applied to global galaxy populations (eg., cosmic star 
formation history by Madau, Pozzeti & Dickinson 1998), 
but this is the first application specifically targeting cluster 
galaxies. We also address how the blue fraction varies with 
galaxy magnitude and with radius from the cluster centre. 
The variation with magnitude is of interest given Cowie's 
claim of cosmic 'down-sizing': the progression of star for- 
mation into smaller and smaller units as the universe ages 
(Cowie et al. 1996). 

Clearly, the history of galaxy accretion will vary from 
cluster to cluster, and this will limit the degree to which we 
can reconstruct cluster galaxy evolution from a set of obser- 
vations of different clusters. However, it is important to dis- 
tinguish between the growth of the total cluster mass, which 
will be dominated by mergers between similar mass units, 
and the accretion of field galaxies. Monte-Carlo simulation 
of cluster growth (Lacey & Cole, 1993) shows that while the 
sequence of mergers varies greatly between clusters, the role 
of accretion is more stable, and the fraction of galaxies ac- 
creted varies ~ 2 per cent per Gyr. This is particularly the 
case if we base our study on X-ray selected clusters since the 
existence of a dense intracluster medium biases against the 
inclusion of clusters which are not yet dynamically relaxed 
(Bower et al. 1997). Nevertheless, differences in cluster ac- 
cretion histories can only be eliminated using a much larger 
sample of clusters than we have available at present. 

It is also important that we try to match the masses 
of the clusters in order to form an evolutionary sequence. 
The X-ray luminosities of the CNOC clusters used in this 
study bracket that of the Coma cluster (§ 2.1). Since the 
X-ray luminosity function evolves little over the look-back 
times involved (eg., Ebelling et al. 1996), the clusters should 
be reasonably well matched in terms of their space density. 
Possibly, we would have liked to compare with an even richer 
local cluster (§ 2.4). However, although there is no such mas- 
sive local cluster, current data on the richness dependence of 
the CM relation suggests that it would have even an tighter 
CM ridge-line (and smaller blue fraction) than that of Coma 
(Lopez-Cruz 1996). 

Our goals in this paper are two-fold. Initially, we will 
make an accurate comparison of the distant and local clus- 
ters. In particular, we will test whether the galaxies in the 
distant clusters simply fade and become red after their star 
formation is truncated or whether an additional mechanism 
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such as tidal stripping or harassment (Moore et al. 1996) is 
required to explain how the distant clusters evolve into clus- 
ters like those at the present-day. We then go on to explore 
the distribution of galaxies in the colour-magnitude plane, 
and to address the difference between clusters at different 
redshifts. The distribution of galaxies can be used to derive 
the history of galaxy accretion from the field and the way the 
star formation is truncated. We will show that a consistent 
picture emerges in which the accretion rate of cluster galax- 
ies from the field declines smoothly with time. Combined 
with the decline of star formation activity of field galaxies, 
this naturally accounts for the Butcher-Oemler effect. 

The layout of the paper is as follows. In § 2, we present 
the sources of data used to compose the distant and local 
cluster CM diagrams, and outline how we correct for back- 
ground and foreground galaxies. In § 3, we calculate the blue 
galaxy fraction from the field corrected CM diagrams and 
discuss its dependence on magnitude and the distance from 
cluster centre. In § 4, we apply our stellar population models 
to calculate the colour change (red-wards) and fading of the 
distant cluster CM diagrams and compare with a local clus- 
ter (Coma). In § 5, we investigate the colour distribution of 
galaxies in more detail, examining how the galaxy accretion 
rate must evolve in time. Our conclusions are summarised 
in § 6. 

We use Ho=50 km s _1 Mpc -1 and qo=0.1 throughout 
this paper, which are the same parameters adopted in the 
original B084 work. The age of the universe is 16.5 Gyr for 
this parameter set. 



2 PHOTOMETRIC DATA FOR CLUSTERS OF 
GALAXIES 

2.1 Distant Clusters 

We have chosen to use the CNOC clusters (Yee et al. 1996), 
since they form a homogeneous set of intermediate redshift 
clusters, selected, with one exception, from the EMSS cata- 
logue of X-ray bright clusters (Gioia & Luppino 1994). Al- 
though the CNOC clusters consist of 16 clusters in total 
covering 0.17<2<0.55, we have restricted attention to the 
redshift range 0.23^z<0.43. We have set the lower limit for 
the cluster redshift so that the CNOC photometry data in 
Gunn g and r (effective wavelengths 4900 A and 6500 A, re- 
spectively) bracket the 4000-A break. Small and well-defined 
k-corrections then allow us to reliably compare the clusters 
in rest frame U—V colours. We set an upper-limit for the 
cluster redshift so that the CNOC photometry is accurate 
enough to recover reliable information on galaxy stellar pop- 
ulations. At z ~ 0.55, for example, the photometric error is 
as high as 0.2 magnitude even for the relatively bright clus- 
ter members (~M,). There are 9 clusters within these red- 
shift limits, but the photometric catalogues are not available 
for two clusters (MS1231, MS1455) at z ~ 0.25. The final 
subset of 7 clusters is summarised in Table 1. The clusters 
are all X-ray bright with L x = 4 - 28 x 10 44 erg s" 1 (0.3- 
3.5 keV), which is comparable to the nearby Coma cluster 
(L x = 8.1 x 10 44 erg s" 1 ) (Yee et al. 1996; Ebeling et al. 
1996; see Table 1 of this paper). The advantage of using the 
CNOC sample is that they obtained ~100-300 redshifts per 
cluster and also the field of view of both the imaging and the 



spectroscopy is relatively large, particularly for four clusters 
(Table 1). Both of these help us in determining cluster mem- 
bership accurately (see Appendix A). Furthermore, the wide 
field enables us to investigate systematic changes of stellar 
populations as a function of distance from the cluster centre 
(eg., Abraham et al. 1996). 

We use the positions, colours and magnitudes in g and 
r bands, and spectroscopic redshifts of galaxies, taken from 
the photometry catalogs provided by the CNOC group (Yee 
et al., 1996). The photometry catalog uses a diameter aper- 
ture of £6", which corresponds to £29-43 kpc for 2=0.23- 
0.43. The completeness limit is about r=23.5 magnitude 
(Yee & Ellingson 1996), and we use only those galaxies 
brighter than this limit in this paper. Even so, the pho- 
tometric data has relatively poor accuracy compared to the 
data available for local clusters. The typical photometric er- 
ror reaches to ~0.2 magnitude at r~23 mag. 

2.2 Matching the Cluster Radius 

One of the key ingredients required to make the accurate 
comparison of galaxies between local and distant clusters is 
to match the area over which the comparison is made. We 
use i?3o , the radius from the cluster centre which contains 30 
per cent of the whole cluster populations, following B084. 

R30 is determined using the galaxy distribution from 
the cluster centre after subtracting the average number den- 
sity of field galaxies taken from the edge of the cluster field 
(see below). Since the shape of the observed field is rectan- 
gular, we make up the galaxy distribution outside of the 
observed region by re-sampling the galaxies at the same 
cluster-centric distance. We measured R30 for four clusters, 
namely, A2390, MS1358+62, MS1512+36, and MS1621+26. 
These are listed in Table 1, together with the concentration 
index C: 

C = log(J?6o/ii2o), (1) 

where R20 and 7?60 are the radii which contain 20 and 60 
per cent of the whole cluster populations, respectively. C 
is always larger than 0.35, and they are classified as either 
compact or intermediate clusters. Our estimates are in good 
agreement with those of van der Marel et al. (1999) ex- 
cept for MS1512+36. As we will discuss later, MS1512+36 is 
likely to be contaminated by a foreground group or cluster, 
making the -R30 estimate difficult. For the remaining three 
clusters, it is impossible to obtain R30 because of their lim- 
ited spatial coverage. Therefore we assume R30 to be 1 Mpc 
(Ho =50), the typical radius that we measured for the other 
clusters. This may be in error by ~30 per cent due to fluctu- 
ation of the background field galaxy density. We will treat 
these clusters separately in our discussion. Note that -R30 
always falls within the original field of view. 

2.3 Field Galaxy Subtraction 

The next ingredient that is needed for the accurate com- 
parison between local and distant clusters is the fore- 
ground/background galaxy subtraction. We use the spectro- 
scopic redshifts as the primary information cluster member- 
ship, where they are available. We define the cluster mem- 
bers as those whose peculiar velocities are within ± 3000 km 
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Table 1. A summary of the cluster sample. Each row shows (1) the cluster name, (2) redshift, (3) X-ray luminosity in the 0.3-3.5keV 
band, (4) field coverage, (5) look back time at cluster redshift, (6) radial cut R30, (7) B084 magnitude cut (My=— 20), (8) slope of the 
CM sequence, (9) zero-point of the CM sequence at r cu t — 1.5, (10) B084 criterion for the blue galaxies defined as the distance from the 
CM sequence (A(B— V)=— 0.2), and (11) fraction of blue galaxies and Poisson error (1<t). 
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A2390 


0.228 


23.8 


46 x 7 


3.59 


2.77 


20.90 


-0.020 


0.89 


-0.26 


0.105 ± 0.034 


MS1008-12 


0.306 


4.5 


9.0 x 7.9 


4.52 


2.82+ 


21.71 


-0.029 


1.20 


-0.36 


0.161 ± 0.036 


MS1224+20 


0.320 


4.6 


9.0 x 7.1 


4.67 


2.74t 


21.84 


-0.031 


1.12 


-0.37 


0.174 ± 0.058 


MS1358+62 


0.327 


10.6 


9.0 x 23.3 


4.75 


3.58 


21.90 


-0.032 


1.16 


-0.38 


0.106 ± 0.029 


MS1512+36 


0.371 


4.8 


27.3 x 8.0 


5.21 


3.26 


22.27 


-0.037 


1.18 


-0.40 


0.379 ± 0.074 


MS0302+17 


0.425 


5.0 


9.0 x 7.9 


5.73 


2.31+ 


22.72 


-0.041 


1.46 


-0.41 


0.213 ± 0.075 


MS1621+26 


0.428 


4.5 


9.0 x 23.3 


5.75 


2.18 


22.75 


-0.041 


1.64 


-0.40 


0.159 ± 0.039 


Coma 


0.024 


8.1 


63 x 53 


0.46 


22.0 










0.056 ± 0.028 



Note - t R 30 is assumed to 1 Mpc for these clusters due to the relatively small spatial coverage. 



s _1 from the cluster central velocity. Galaxies outside of this 
redshift range are assigned to the 'field'. However, the depth 
of the spectroscopic redshifts is quite shallow (r=21.5), and 
there are many galaxies for which spectroscopic redshifts are 
unavailable even brighter than this magnitude limit. 

We use a statistical method to determine cluster mem- 
bership for the galaxies without spectroscopic information. 
The details of this method are given in Appendix A and 
we give a brief summary below. We take the sample of field 
galaxies from the edge of the clusters A2390, MS1358+62 
and MS1621+26. From the distribution of these galaxies on 
the CM diagram, we determine the probability that a given 
galaxy in a cluster with a given colour and magnitude is a 
field galaxy. Using this probability, we apply Monte-Carlo 
simulations to assign cluster membership statistically and 
to subtract the statistical non-members from the CM dia- 
grams. The Monte-Carlo simulation is repeated 100 times 
per cluster and we average the realisations for the analy- 
ses presented in this paper. Filled symbols in Figs. 1 and 2 
show the resulting field corrected CM diagrams for a typical 
realisation. 



2.4 Local Cluster Comparison 

We use Coma as a representative rich local cluster for com- 
parison with the CNOC clusters at higher redshifts. We take 
the high precision Coma photometric data from Terlevich 
ct al. (2000). They imaged Coma in U and V bands with 
a spatial coverage of 3360 arcmin 2 . We use a fixed diam- 
eter aperture of 25.3 arcsec, both for the magnitudes and 
the colours, which corresponds to a physical scale of 17 kpc. 
The measurement error is very small, and the typical error 
in U—V colour is less than 0.035 magnitude. R30 for Coma is 
22 arcmin or 0.9 Mpc (B084; Merritt 1987). We fit the CM 
relation for galaxies within R30 and brighter than Mv>— 20 
using the Bi-weight method. We also define a 'red envelope' 
at 0.073 magnitude redder than the fitted CM relation (red 
sequence). This colour offset corresponds to a colour differ- 
ence of the SSP models between 2f orm =2 and 5.4 with solar 
metallicity. The CM diagram of Coma within R30 is given 
in the bottom right panel of Fig. 7. 

Note that we will not take into account the dynamical 
evolution of clusters in this paper. An estimate of this effect 



shows that it is unlikely to be important if we restrict at- 
tention to the region within R30 ■ According to the extended 
Press-Schechter (P-S) theory (Bower 1991; Bond et al. 1991; 
Lacey & Cole 1993), a rich cluster comparable to Coma had 
acquired ~60 per cent of the present-day mass by 2=0.33 
(qo=0.1). Judging from the X-ray luminosity (Table 1), two 
of our sample clusters, A2390 and MS1358+62, will evolve 
into richer systems than Coma before the present-day. For 
these systems, we would have liked to compare with an even 
richer local cluster. Although there is no such rich coun- 
terpart nearby, data on the dependence of the CM relation 
with richness suggests that it would have even an tighter 
CM ridge-line than that of Coma (Bower et al. 1992; Lopez- 
Cruz 1996). If this is the case, the recent star formation 
would be more restricted and earlier truncation would be 
required for these clusters than that discussed in this paper. 
In addition, most of the galaxies within R30 at intermedi- 
ate redshifts will appear in slightly smaller regions than R30 
at the present-day, because the cluster will acquire another 
^40 per cent of new galaxies between 2=0.33 and 2=0 for 
example (Bower 1991). Therefore, the population within R30 
at a 2=0.33 cluster should be actually compared with those 
within Ris at 2=0, where Ri$ is the radius which contains 
18 per cent of the whole cluster population. However, the 
colour scatter within R30 of Coma is already so small that 
it hardly changes within this radius. 



3 THE BUTCHER-OEMLER EFFECT 
3.1 Blue Galaxy Fraction 

From the field corrected CM diagrams (Figs. 1 and 2) we now 
calculate the blue galaxy fraction in each cluster. We follow 
the original B084 definition of the blue galaxy fraction: only 
galaxies within 7?3o and brighter than My- =—20 in the rest 
frame are considered; only those bluer than A(B — V)=— 0.2 
from the red sequence are defined as the 'blue' cluster mem- 
bers. We transform these criteria for blue galaxies into the 
ones for the CNOC g and r bands in the observer's frame 
using the K-corrections given in Fukugita, Shimasaku & 
Ichikawa (1995). They are based on Kennicutt's (1992a) 
and Coleman, Wu & Weedman's (1980) spectral atlas. Ac- 
cording to Fukugita et al. (1995), A(B-V)=-0.2 roughly 
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Figure 1. Colour-magnitude diagrams for A2390 (top left), MS1008-12 (top right), MS1224+20 (bottom left), and MS1358+62 (bottom 
right) within -R30. The symbols are as follows:- filled circles, spectroscopically confirmed members; crosses, spcctroscopically confirmed 
non-members; filled triangles, statistically plausible members; open triangles, statistically subtracted non-members. The lines in the 
figure show:- horizontal long dashed line, the predicted CM red sequence; solid line, the red envelope (see text for definition); dotted line, 
colour limit of the blue galaxies corresponding to B084's limit of A(£?— V)=— 0.2 from the red sequence; vertical dashed line, magnitude 
cut, r cu t, corresponding to My=— 20 in the rest frame as in B084. Error bars in the figure illustrate typical la photometric errors and 
their dependence on magnitude. 



corresponds to the colour of Sab-type galaxies in the Hub- 
ble sequence of the present-day galaxies, assuming that the 
zero-point corresponds to the E-type colour, and it can be 
transformed into A(g— r) at a given redshift using the colour 
difference between E and Sab types. The B084 criteria 
A(B — V)=— 0.2 actually falls between Sab and Sbc types, 



and we make an interpolation to get the exact A(g— r). This 
also corresponds to A(U—V)=0.54 at the present-day, which 
is used as a B084 criterion for our Coma data. The mag- 
nitude cut, Mv=—2Q, is also transformed into an apparent 
magnitude in r-band, r cut , for each cluster, assuming the 
E-type spectral energy distribution (SED) as in B084. The 
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Figure 2. Colour-magnitude diagrams for MS1512+36 (top left), MS0302+17 (top right), and MS1621+26 (bottom left) within R 30 , 
and for the field taken from the edge of the three cluster fields (bottom right; see text). The symbols and lines are explained in the 
caption of Fig. 1. 



adopted r cut and A(g—r) for each cluster are summarised in 
Table 1. These limits are also shown by the dot-dashed and 
the dotted lines, respectively, in Figs. 1 and 2. 

However, the above magnitude cut is not quite correct 
for the blue galaxies. A blue galaxy with a given observed r 
magnitude is actually fainter in Mv than a red galaxy with 
the same r magnitude if the r-band in the observed frame 
corresponds to shorter wavelength than the rest frame V- 
band, as is actually the case for our CNOC clusters. There- 



fore, we apply colour dependent K-corrections using the 
numbers for various spectral types given in Fukugita et al. 
(1995). The B084 magnitude cut in r-band thus determined 
is shown by the slanted dashed line in Figs. 1 and 2. The 
deviation from the constant r cu t (dot-dashed) is larger for 
blue g— r colours and increases with redshift. The simpler 
definition of the magnitude cut (as is adopted in BO) tends 
to overestimate the fraction of blue galaxies by a few per 
cent (A/ s ~0.2) (see Fig.3). 
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For each cluster, we define the location of the red se- 
quence on the CM diagram as a zero-point from which we 
measure the colour difference. Since the CM slope is poorly 
determined with the CNOC photometry, we use a model 
prediction, assuming the CM slope is governed by metallic- 
ity variation. We adopted the metallicity sequence model in 
Kodama & Arimoto (1997). The model slope is calibrated to 
reproduce the Coma CM slopes in U—V and V—K (Bower 
ct al. 1992) at 15 Gyr. Since the CM slope depends on the 
aperture within which we integrate colours, we adopt the 
50 kpc aperture models to closely match the CNOC's large 
aperture (cf., Kodama et al. 1998). Using this fixed slope, 
the CM red sequence is fitted to the observed data points so 
that the dispersion around the fit is minimised after remov- 
ing the galaxies that are located outside of the one sigma 
limits for all the galaxies. The slopes and the zero-points at 
'/"cut — 1-5 (Mv=— 21.5) of the adopted CM red sequences are 
given in Table 1. 

We also define a 'red envelope' which will be used later 
in § 4. We assume that the 'red sequence' corresponds to 
Zform=2, supposing it is the typical star formation epoch of 
cluster ellipticals (eg., Bower et al. 1992; Ellis et al. 1997; 
Stanford, Eisenhardt & Dickinson 1998; BKT98; Kodama 
et al. 1998; van Dokkum et al. 1998a; 1998b), and draw 
the red envelope parallel to the red sequence so that the 
envelope is redder than the sequence by the same amount of 
colour difference (in g—r) between a SSP model of Zf orm =2 
(T G =12.3 Gyr) and that of z form =5.4 (T G =15 Gyr). The 
colour difference ranges from 0.05-0.1 magnitude. 

Following all these steps, we can now reproduce the def- 
inition of the blue galaxy fraction (/s) used by B084, with 
the addition of a colour dependent magnitude cut as dis- 
cussed above. We determine error bars by combining the un- 
certainty in the blue fraction from the statistical subtraction 
with a Poisson error based on the numbers of blue galaxies: 



total 



^cluster + N] 



field, stat 
blue 



^.cluster 
Afield 



(2) 



where iVtotai is the total number of galaxies after field 
subtraction, A^f 01 is the number of blue members, and 
^biuo ' Stat ^ s the number of blue field galaxies that are sub- 
tracted statistically, A c i us t cr and ^flcw are the area of the 
cluster and the area of the field sample, respectively. The 
Poisson component always dominates the error. The blue 
fraction of each of the clusters is summarised in Table 1. 

For A2390 and MS1621+26, we can cross-compare $b 
with the earlier estimates by Abraham et al. (1996) and 
Morris et al. (1998), respectively. / B =0.14±0.05 for A2390 
and /s~0.2 for MS1621+26 at our R30 are both consistent 
with our estimate within the errors. We can also compare our 
results for A2390 with the study of a larger sample of X-ray 
luminous (L x (0. 1-2.4) > 16 • 10 44 ergs" 1 ) clusters at z~0.2 
by Smail et al. (1998). Working wih rest-frame B— I colours, 
they found an average blue fraction of 0.032 ± 0.069. Our 
estimate for A2390 is somewhat larger than this number, 
but consistent at the 1 a level. 

Figure 3 shows /s as a function of cluster redshift. Our 
results fit the original evolution line (solid line) suggested 
by B084 remarkably well, ft is notable, however, that one 
cluster, MS1512+36 at 2=0.37, has anomalously high blue 
galaxy fraction compared to the B084 line and the other 
CNOC clusters. The CM diagram of this cluster suggests 




0.2 0.3 0.4 

Redshift (z) 



0.6 



Figure 3. The blue galaxy fraction (3b) in clusters as a func- 
tion of redshift (see text for definition). The error bars indicate 
the Poisson errors (1<t). The solid line shows the evolution line 
originally suggested by B084. Open circles show the clusters for 
which we had to assume R30 due to the small spatial coverage. 
The Coma cluster is plotted as a filled square. Crosses show the 
effect of neglecting the colour dependent K-corrections for the 
magnitude cut (see text for detail). MS1512+36 in the open tri- 
angle shows an exceptionally high fraction of blue galaxies, but 
it is likely to be contaminated by a superposition of a foreground 
group/cluster (see text). 



that there may be a foreground group/cluster superimposed 
along the line of sight. As seen in the top left panel of 
Fig. 2, there are ~5 bright galaxies (r<19) with similar 
colour (g— r~0.6) within R30 which are not subtracted in 
the statistical realisations. This clump appears to continue 
down to fainter magnitudes and to form another CM se- 
quence. Although most of the galaxies on this secondary 
sequence do not have spectroscopic redshifts (which makes 
it difficult to identify the foreground structure conclusively) 
there are 5 more galaxies in the entire cluster field which 
have similar colours of g— r~0.6 and spectroscopic redshifts 
in the narrow range 0.163T0.002. At this redshift, the colour 
is consistent with passively evolving early-type galaxies (Ko- 
dama & Arimoto 1997), which leads us to suggest that the 
high value of the blue fraction of MS1512+36 is contami- 
nated by the foreground group/cluster. 



3.2 Dependence on Radius 

One of the interesting questions concerning the blue galaxies 
is their spatial distribution. If the galaxies are continuously 
accreted from the surrounding field, we should expect the 
blue galaxies to be located in the outskirts of the cluster. 
To test for this radial dependence, we separate the CM di- 
agrams into 4 bins within 2 i?3o. If the radius exceeds the 
shorter side of the rectangular observed strip, we make up 
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Figure 4. Dependence of the blue galaxy fraction on the radius 
(R) from the cluster centre. The radius is normalised using R3Q. 
The same symbols are used as in Fig. 3 depending on the cluster. 
The morphology-radius (T-R) relations at z~0 and ~0.5 (Dressier 
et al. 1997) are plotted in the bottom middle and the bottom 
right panel, respectively, assuming that the spiral galaxy fraction 
is directly connected to the blue fraction and that i?3o=l Mpc. 

the galaxy population outside of the strip by statistically 
sampling the galaxies at the same radii within the observed 
field. 

Fig. 4 shows a general trend for f B to increase with 
radius. This is consistent with the previous results shown 
by B084, Abraham et al. (1996) and Morris et al. (1998). 
This trend is to be expected, since it mainly reflects 
the morphology- density (or morphology-radius) relation 
(Dressier et al. 1997) which shows an increasing (decreas- 
ing) fraction of spiral (elliptical) galaxies with galaxy density 
or radius. For comparison, we reproduce the morphology- 
radius (T-R) relation by Dressier (1997) for the local clus- 
ters (z~0) and for the intermediate redshift clusters (2~0.5) 
in the bottom middle and the bottom right panels of Fig. 4, 
respectively. Our sample clusters fall between these two red- 
shifts. Here we assume the blue fraction is equivalent to the 
spiral fraction, since the B084 colour cut for the blue galax- 
ies corresponds to Sab type in the local universe as already 
mentioned earlier (§ 3.1; Fukugita et al. 1995). The slope of 
these T-R relations are steep enough to explain the radial 
dependence of the BO blue fraction. This radial dependence 
also illustrates how important it is that distant and local 
clusters are compared within equivalent radii. 

3.3 Dependence on Magnitude 

Another important issue is the magnitude dependence of 
the blue fraction as shown in Fig. 5. In most of the clus- 
ters, there is only a weak trend for the blue galaxy fraction 
to increase towards fainter magnitudes. MS1512+36, which 



Figure 5. The dependence of the blue galaxy fraction on abso- 
lute magnitude. The fiducial magnitude cut for fi corresponds to 
My=— 20. The same symbols are described in Fig. 3. The points 
at the faintest magnitude bins for the two highest redshift clus- 
ters are dominated by the large photometric errors (~0.2) and 
are not shown. 

has the strongest trend, is likely to be contaminated by a 
foreground group/cluster (§ 3.1). The overall impression is 
that the magnitude dependence is weak especially for the 
bright magnitudes (My>— 20). B084 found no systematic 
change in blue fraction over a similar magnitude range, ex- 
cept for their brightest magnitude bin (— 23<My<— 22.5), 
where the blue fraction is smaller. However, these results do 
not mean that there is little mass dependence of the star for- 
mation histories (Smail et al., 1998, Poggianti et al., 1999). 
In order to discuss the intrinsic mass dependence we need to 
consider the fading of the blue galaxies after they enter the 
cluster, since the star forming galaxies are brighter in My 
than the red galaxies of the same mass. We will come back 
to this point later in § 4.3 by showing the properties of the 
galaxies as a function of the 'present-day' magnitude. This 
approach will give us a direct view to the mass dependence 
of star formation history. 

4 FROM INTERMEDIATE REDSHIFT 
CLUSTERS TO THE COMA CLUSTER 

In this section, we address the fate of the blue galaxies seen 
in the intermediate redshift clusters. In particular, we inves- 
tigate whether they become incorporated into the tight CM 
relation seen in the present-day local clusters. 

BKT98 discussed this issue, introducing the models 
where the star formation is truncated at the time of galaxy 
accretion from the surrounding field to clusters. By estimat- 
ing the distribution of the truncation epoch by the extended 
P-S theory, they explored the allowed range of star forma- 
tion histories in rich clusters. In this paper, we discuss this 
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problem in further detail, taking into account the fading of 
the galaxies and their colour evolution after star formation 
ceases. Our approach extends the analysis of Smail et al. 
(1998) who addressed the evolution of galaxies in a sam- 
ple of X-ray luminous clusters at z ~ 0.2. Instead of using 
the predicted galaxy accretion rate at the first place, we 
start from the real galaxy distribution on the CM diagrams 
of intermediate redshift clusters, apply the star formation 
truncation model for the individual galaxies, and sketch the 
whole movement of cluster galaxies on the CM diagrams 
down to the present-day. Our approach does not assume a 
particular galaxy accretion history for the clusters: we dis- 
cuss the limits that can be placed on the accretion history 
in §5. 

In the following sections, we concentrate on the results 
for the three clusters for which the background subtraction 
is most reliable: A2390, MS1358+62, and MS1621+26, (see 
§ 2.3 and Appendix A). MS1512+36 has been omitted be- 
cause of the probable contamination by a foreground sys- 
tem. The three other clusters from the CNOC sample have 
only a relatively small area of imaging making the back- 
ground subtraction unreliable. Nevertheless the results for 
these clusters agree qualitatively with the systems on which 
we now concentrate. 

4.1 Modelling Galaxy Evolution Vectors 

Initially, a difficulty with this approach appears to be that 
the photometric evolution of a galaxy after the end of its 
star formation, will depend on its star formation history 
prior to the event. Fortunately, this is not the case: so long 
as a fraction of stars are old, the passive evolution phase 
is only weakly dependent on the previous history. Fig. 6 
shows some examples of evolution vectors from 2=0.33 to 
for three galaxies with different colours (g— r=0.6, 0.9, 
and 1.2). For each galaxy, we plot 5 different vectors which 
correspond to different star formation histories. Star forma- 
tion is supressed prior to z=0.33 for all the models, with the 
truncation time (ttrunc) being chosen so that each model has 
the same g— r colour at the moment of observation. Apart 
from the single burst model which gives the maximum fad- 
ing and the bluest colour at the present-day, all the other 
models show similar evolution vectors on the CM diagram 
in spite of large variations in their previous star formation 
histories. This remains true even if a burst of star formation 
is included just before the truncation epoch. As a result, 
the evolution of the cluster members on the CM diagrams 
depends little on their prior star formation history. 

We use the gas infall model (described in detail in Ap- 
pendix B) to parameterise the star formation of a galaxy 
before it enters the cluster environment. In outline, gas in- 
flow onto the galaxy determines its star formation timescale 
and hence the distribution of field galaxy colours. The dis- 
tribution of gas infall timescales is set so as to match the 
observed colours of the field galaxies in the CNOC sam- 
ple and the global star formation in the Universe (Madau 
ct al. 1998). Once a galaxy is accreted by the cluster, star 
formation is truncated, possibly after a short burst (as de- 
scribed in Appendix C). Adding the burst helps to explain 
the frequent appearance of strong Balmer absorption fea- 
tures in the integrated spectra of galaxies in the intermediate 
redshift clusters (Barger et al. 1994; Poggianti et al. 1999; 
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Figure 6. Evolution vectors, showing the fading and colour evo- 
lution after star formation ceases, in the rest frame U—V vs. 
My diagram. The length of the vector shows the evolution ex- 
pected between 2=0.33 and z=0. Three sets of vectors are plotted 
to show the evolution of three galaxies with different colours at 
z=0.33; (r,g-r)=(21,0.6), (21,0.9), and (21,1.2). For each galaxy 
five vectors are shown, correspond to different star formation his- 
tories prior to the truncation event:- (1) dotted line, constant star 
formation; (2) short dashed line, exponentially decaying star for- 
mation with an c-folding time of r=5 Gyr (simple model); (3) 
solid line: gas infall model with t=3 Gyr (see Appendix B for 
details); (4) long dashed line, as (3), but with a 10 per cent star 
burst added just before truncation; (5) dash-dotted line, single 
star burst just before the truncation. We assume that the star 
formation starts at £lb=15 Gyr (z=5.4) for all the models except 
(5) and terminates abruptly. If star formation decays slowly after 
the truncation event, the direction and the length of the vectors 
is hardly changed, although the speed of evolution along the vec- 
tor immediately after truncation does change significantly. The 
horizontal solid line shows the red envelope at z=0. 



BKT98; but see also Balogh et al. 1999). Although there 
is ambiguity in the star formation and the gas infall time 
scales and the burst strength here, the subsequent evolution 
is rather insensitive to these values unless the burst popula- 
tion dominates. We assign the star formation/gas infall time 
scale r and the truncation time ttrunc for each galaxy to re- 
produce its colour in g—r at the cluster redshift as explained 
in Appendix D in detail. Our approach differs from that of 
Smail et al. (1998) who applied a similar kind of trunca- 
tion model, under the assumption that star formation was 
terminated in all the cluster members at the same epoch. 
Instead, we allow ttrunc to vary from galaxy to galaxy re- 
flecting the way the cluster is being built up continuously 
over time. With this parameterisation, we can later discuss 
the galaxy accretion history which forms the clusters (§ 5). 
Once the r and ttrunc is assigned for each cluster member, 
the model gives its colour and magnitude evolution in rest 
frame U— V colours and absolute My magnitudes. 
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The average metallicity of stellar populations varies 
with magnitude along the CM sequence (Kodama et al. 
1998; Terlevich et al. 1999; Kuntschner 2000). This varia- 
tion can either be incorporated into the model fading vec- 
tors, or the observed galaxy colours can be corrected to a 
fiducial metal abundance (as in BKT98). Since the fading 
vectors depend little on metallicity around the solar value 
(Z=0.02), it is simpler to adopt the second approach. There- 
fore, we corrected the observed colours for the CM slope as 
follows: 

{g - T - ) corrected = (3 ~ f") observed ^ sl °P C X i T ~ r °) ( 3 ) 

where ro=r cu t — 1-5 which corresponds to My=— 21.5 in the 
rest frame. Note that although we can correct the systematic 
trend as a function of luminosity, we neglect the effect of the 
age- metallicity conspiracy (Worthey, Trager & Faber 1996; 
Trager 1997; Shioya & Bekki 1998) in which young galaxies 
are supposed to be more metal rich (and hence the age and 
metallicity effects on galaxy colours may cancel out each 
other to some extent). Even if this is the case, this effect 
would be important only for the red galaxies on the CM 
sequence. We will return to consider the implications of the 
CNOC clusters for age variations along the CM relation in 
§ 4.3. 

Before applying the truncation model, we have to cor- 
rect a small zero-point mismatch between the observed 
colours and the model by comparing between the colour of 
the red sequence at ro and the SSP model with solar metal- 
licity and Zform=2. We shifted the observed g— r colours (al- 
ready corrected for the CM slope) by 0.0-0.15 magnitude. 
Although this is probably due to the zero-point uncertainties 
in the model and the observation, it may also be due to the 
metallicity mismatch since we fixed the model metallicity at 
solar metallicity. 

A further problem is that there are a few galaxies which 
have been statistically assigned cluster membership, even 
though they are redder than the CM sequence. These objects 
cause problems because they are difficult to model within 
our stellar population code. Some of these galaxies may be 
background objects which have not been removed by the 
statistical subtraction process. Alternatively, these red ob- 
jects are likely to be scattered from the red sequence due to 
the relatively large photometric errors (~0.2 mag) toward 
fainter magnitudes. Unlike galaxies scattered blueward of 
the sequence, the offset of these galaxies would remain large 
as we evolve the cluster forward in time, yet there are no 
such galaxies seen in the Coma cluster. We have therefore 
adopted the scheme of initially correcting these red galaxies 
back onto the CM relation sequence. A similar correction 
is not required for the blue galaxies since the differences 
due to photometric errors decline as the stellar population 
ages. It should be noted that the scatter of the red ob- 
jects from the CM sequence is not large enough to affect 
the blue galaxy fraction (/s) significantly — as can be seen 
by comparing the deviation of red and blue objects from 
the CM sequence in Figs. 1 and 2. Most of the red objects 
fall within A(g— r)<0.4 from the CM red envelope, of which 
most of the spectroscopically confirmed members are within 
A(gr— r)<0.2. It is possible these red galaxies are dusty spi- 
rals. If these have extensive on-going star formation, we will 
tend to underestimate the fraction of star forming galaxies 



in the cluster. We discuss the consequences of dusty objects 
within clusters in § 6. 

4.2 Evolving the Distant Clusters Forward in 
Time 

We are now in a position to simulate the photometric evo- 
lution of the cluster galaxy population. Fig. 7 illustrates the 
evolution of the three clusters forward in time down to the 
present-day (from the left panel to the right). All the cluster 
members within R30 are plotted after the field subtraction. 

Figure 7 is an example of the Monte-Carlo simulations 
using ei oss =0 and eb urs t=0.01 (Model (b) in Appendix C). 
The colour is presented with respect to the red envelope at 
2=0 as discussed in §4.1. Note that as we evolve the cluster 
forward in time (the panels to the right of the initial data) , 
there is no further accretion. Within a vertical row, the dif- 
ferences between the panels reflect the accretion that has 
occured in the younger cluster. 

If we compare the initial CM diagrams (left end panels) 
of those intermediate redshift clusters with that of the Coma 
cluster (bottom right panel) in the same coordinates, it is 
obvious that there are more blue galaxies in the distant clus- 
ters and the scatter around the CM red sequence is larger. 
However, if we trace the movement of those blue galaxies for- 
ward in time, we see that the galaxies become redder rapidly 
as time passes and all get redder than the B084 criteria of 
blue galaxies (dotted line) in only a few Gyr at most after 
the cluster redshifts. As a result, the colour scatter of the 
CM relation becomes small enough by the present-day, es- 
pecially at the bright end (Mv<— 20), to be consistent with 
the observed scatter in Coma. A more precise comparison is 
not justified by the photometric accuracy of the CNOC data 
and our consequent treatment of objects scattered redward 
of the initial CM relation. If galaxy harrassment accompa- 
nies the galaxy accretion, stripping off some stars from the 
galaxy (as suggested by numerical simulations: Moore et al. 
1996; Moore et al. 1999a), the blue galaxies could fade even 
further, resulting in a smaller present-day colour scatter at 
the bright end of the CM sequence. On the basis of this com- 
parison, such an effect is not required by the data, however. 

4.3 The Mass Dependence of Star Formation 
Histories 

Not only do the galaxies become redder, but they also get 
fainter after the truncation event. In fact, most of the blue 
galaxies which pass the B084 criterion at the cluster red- 
shifts fade rapidly and evolve to red faint galaxies with 
Mv>— 20 (M,+l) by the present-day as shown in Fig. 8 
(Fig. 7 shows the trajectory of galaxies on the CM dia- 
grams). The fraction of blue galaxies is now plotted as a 
function of the present-day magnitude instead of the mag- 
nitude at the observed epoch (as in Fig. 5). Therefore, the 
blue galaxies at intermediate redshifts are not the counter- 
parts of the bright early-type galaxies seen in present-day 
clusters including SO's (Smail et al., 1998). This is consis- 
tent with Kelson et al.'s (1999) and Poggianti et al.'s (1999) 
result derived from the spectroscopic signatures of star for- 
mation activity. It would appear that bright (hence massive) 
early-type galaxies should have completed their star forma- 
tion well beyond 2=0.4, whereas fainter (hence less massive) 
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Figure 7. Evolution of the cluster galaxies within R30 forward in time on the CM diagrams. From the left panel to the right, redshift 
is decreasing as 0.43, 0.33, 0.23 and (an evolutionary sequence). Star formation is truncated after the accretion of galaxies from the 
field. The truncation epoch (ttrunc) is assigned for each galaxy to reproduce its g—r colour at cluster redshift. This is an example of 
the Monte-Carlo simulations using ei oss =0 and etj Urs t=0.01 (Model (b)). As we propagate the cluster forward in time, no further galaxy 
accretion occurs. The effect of accretion can be seen by comparing the panels vertically (ie., at the same redshift). The Coma cluster is 
plotted in the bottom right panel for comparison with the z=0 models of the CNOC clusters. Colours are presented with respect to the 
red envelope (solid line) to eliminate the metallicity effect. The dotted line and the dashed line represent the B084 criterion of the blue 
galaxies which corresponds to A(B— V)=— 0.2 and My=— 20, respectively. 



galaxies have had more extended star formation. This seems 
to imply that the cosmic down-sizing hypothesis (Cowie et 
al., 1996; Cowie, Songaila & Barger 1999), where star for- 
mation propagates to smaller units as the Universe ages, 
applies to cluster galaxies as well as galaxies in the field. 



By extending this analysis to include higher redshift clus- 
ters, it should become possible to more strongly constrain 
the formation of massive cluster members. 

Since the fraction of recently accreted galaxies increases 
at faint magnitudes, the resulting age bias may contribute to 
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Figure 8. Fraction of galaxies that are identified as blue galax- 
ies at the observed epoch plotted as a function of the present- 
day (z = 0) magnitude. This realisation results from adopting 
Model (b) as described in Appendix C. The points are slightly 
offset horizontally to avoid their overlapping, ie., —0.1 for A2390 
and +0.1 for MS1621+26. 
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Figure 9. Luminosity weighted mean age at the present-day 
as a function of the present-day magnitude. Model (b) is used 
(see § 5.1). The points are slightly offset horizontally to avoid 
their overlapping, ie., -0.15 for Coma, -0.05 for A2390, +0.05 
for MS1358+62, and +0.15 for MS1621+26. Error bars show the 
scatter (1 cr) within each magnitude bin. 



the CM slope. However, if the luminosity (B-band) weighted 
mean age of the galaxies (iiwm) measured at the present-day 
is plotted as a function of galaxy luminosity, the systematic 
trend is very weak; only ~ 1-2 Gyr difference at most over 3- 
4 magnitude range of the present-day CM relation (Fig. 9). 
The lack of magnitude dependence in luminosity weighted 
mean age reflects the fact that ii wm gets old very rapidly 
once the star formation is truncated and is rather insensi- 
tive to a small amount of recent star formation (BKT98). 
However, the systematic age difference is sensitive to the 
CM slope that we have subtracted to correct the metallic- 
ity effect. To see the maximum age variation that can be 
allowed along the CM relation, we repeated the simulation 
without correcting the CM slope. All the systematic effects 
are now attributed to age. Even in this extreme case, the 
age difference along the present-day CM relation is only <, 
3 Gyr, or less than ~20 per cent, and fails to reproduce the 
observed CM slope in Coma: at most only a quarter of the 
aperture corrected CM slope can be accounted for in this 
way. 



5 IMPLICATIONS FOR STAR FORMATION 
HISTORIES IN CLUSTER GALAXIES 

In the previous section, we have modelled the photomet- 
ric evolution of the clusters' galaxy populations in order to 
evolve the system forward in time. In this section, we look 
at the process from the other view-point: using the observed 
colour distribution to infer star formation and truncation 
histories in the cluster environment from high redshifts down 
to the observed epoch. 



5.1 Assigning Star Formation Histories 

We proceed as follows. As a starting point, we adopt a mix 
of galaxy star formation histories that are indicated by ob- 
servations of distant field galaxies: the details are described 
in Appendix D. For each cluster, we are then able to in- 
fer a distribution of truncation times (following the scheme 
outlined in Appendix D) from the observed colour distribu- 
tion. The whole process is treated in a Monte-Carlo sense, 
and 100 realisations are combined to produce the final dis- 
tribution of truncation times. For example, a galaxy with 
only a small blueward deviation from the red CM relation 
sequence might either be the result of a red (low star forma- 
tion rate) field galaxy that arrived in the cluster recently, or 
of a blue (high star formation rate) field galaxy that arrived 
long before the observed epoch. Throughout, we will assume 
that the truncation of star formation is associated with the 
build up of the cluster from accreting field galaxies. Whereas 
we found that the direction of photometric evolution on the 
CM diagrams in the passive evolution phase has little model 
dependence, the time scale of the evolution immediately af- 
ter the truncation epoch (and hence the colour distribution) 
does have considerable dependence on the way the star for- 
mation is truncated (for example the burst strength and the 
level of residual star formation). It is therefore necessary to 
consider different truncation schemes and to compare the 
cluster formation histories of the different clusters. 

We present results based on two particular models. In 
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Model (a), we consider a violent event that might be ex- 
pected to result from ram-pressure stripping (Abadi et al. 
1999; Fujita & Nagashima 1999) or harassment (Moore et 
al. 1996) of the galaxy as it arrives in the cluster. We as- 
sume that half of the gas in the disk is consumed in a strong 
burst of star formation and that the remaining gas is lost 
from the system and no residual star formation is present 
afterwards. This model is intended to reproduce the burst 
models used to explain the spectroscopic properties of E+A 
galaxies (eg., Couch & Sharpies 1987; Poggianti et al. 1999). 
For the second model, Model (b), we assume that there is 
no gas lost from the disk and that only 1 per cent of disk gas 
is converted into stars by a burst. In this case the dominant 
effect of the cluster is to suppress star formation by 'suffoca- 
tion' — removing the halo of material that would otherwise 
have been slowly accreted by the galaxy (Larson, Tinsley & 
Caldwell 1980; Balogh et al. 2000). All the remaining gas 
is consumed exponentially with an e-folding time of 1 Gyr 
(half the average star formation timescale for field galaxies). 

The two models we have chosen are intended to cover 
the range of possible scenarios. In particular, it should be 
noted that if the burst strength is made weaker in Model 
(a), the model struggles to reproduce the population of blue 
galaxies; if the star formation decay timescale is made longer 
in Model (b), it becomes difficult to reproduce a sufficient 
population of red galaxies in the cluster. Nevertheless mod- 
els that lie between the two possibilities selected can suc- 
cessfully reproduce the observed CM diagrams with accre- 
tion histories that are also intermediate between the two 
models. 

5.2 The History of Galaxy Accretion 

We assign star formation histories for individual cluster 
galaxies statistically, based on the observed colour distribu- 
tion, using the field colour distribution to define the initial 
galaxy colours. The resulting distribution of the assigned 
truncation times can now be used to infer the history of 
galaxy accretion from the field into the cluster. Figure 10 
shows the resulting distribution of t tIunz for galaxies brighter 
than My = — 19 at z~0 for the three distant clusters, plus the 
nearby Coma cluster. This magnitude cut ensures that the 
photometry is complete in all of the clusters. The results 
of our Monte-Carlo simulations have been combined into 
three bins (four bins for Coma). Galaxies in which star for- 
mation has been recently truncated can be more accurately 
age-dated than older systems: thus the age ranges of the bins 
increase to reflect the accuracy of the photometric data. The 
two panels differ in the treatment of the truncation event as 
described above. 

As can be seen, the burst model (Model (a)) produces 
a skewed accretion pattern in which the cluster accretion 
rate is always highest at the epoch of observations. This 
is not consistent with the smooth decline in the accretion 
we had expected, nor is it consistent between the different 
clusters. It is unlikely that this peak is caused by an erratic 
accretion history since we see the same pattern in all of the 
clusters. The fundamental problem with Model (a) is that 
the blue galaxies evolve toward the red sequence so rapidly 
that the model struggles to produce the observed numbers of 
blue galaxies. We note, however, in the current model with 
eburst=0.5, the typical burst strength (/burst) at the time of 



truncation is only a few per cent in mass for those galaxies 
just arrived in the clusters. (Even with eb ur st = l, the burst 
strength is typically 5 per cent). In order to avoid the peak 
in the accretion rate at the observed epoch, we could raise 
the burst strength so that the truncation epoch is pushed 
away into the past to get the same colour. However, this 
becomes effective only when we consider very strong bursts 
of at least 20 per cent in mass for all the galaxies. Under the 
current gas infall model scheme (Appendix B), an external 
gas supply is needed to produce such a strong burst, since 
the required amount of gas is not available for most of the 
galaxies. 

The behaviour of Model (a) can be improved if we as- 
sume that galaxies can orbit in the cluster for up to 1 Gyr 
(the cluster crossing time) before undergoing a star burst. 
This tends to smooth out the peak seen in panel (a) because 
we observe some galaxies prior to the truncation event. How- 
ever, the peaks seen just before the observed epoch persist 
since the bin size is already larger (2 Gyr) than the time 
delay. 

In the second panel, we show the decaying model 
(Model (b)). In this model, the history of accretion is much 
more consistent between the different clusters. In all four 
systems the accretion rate declines monotonically with red- 
shift. It is interesting to compare the evolution of the accre- 
tion rate with that predicted on the basis of extended P-S 
theory (Bower, 1991, Bond et al., 1991). As an example, the 
predicted galaxy accretion history for a z=0.33 cluster of 
mass 10 15 Mq is shown by the dashed lines in the figure. 
This has been estimated from the extended P-S theory us- 
ing a threshold halo mass of 5-10 12 M and 5-10 13 M (for 
the right and left peaks respectively) at which the galaxy 
is counted as having made the transition from the field to 
cluster environment (Bower 1991; Kauffmann 1995). As can 
be seen the first of these models fits the data remarkably 
well. It is also notable that the accretion history inferred 
from the Coma cluster matches smoothly onto that inferred 
from the higher redshift systems. This comparison gives us 
fundamental insight into the origin of the Butcher-Oemler 
effect (ie., the differences between high and low redshift clus- 
ters). Both the declining activity of field galaxies, and the 
declining rate at which galaxies are accreted from the 'field' 
environment play an important role in producing the consis- 
tency (both between the accretion histories of the clusters 
and between the reconstructed histories and the theoretical 
model) that we see in this diagram. 

5.3 Comparison with Spectroscopic Studies 

So far, we have concentrated on reproducing a consistent 
accretion pattern between the clusters based on the colour 
distribution of galaxies. In this section, we compare our re- 
sults obtained based on the photometric work with the re- 
cent spectroscopic studies of galaxy evolution in clusters. In 
order to estimate the fraction of the H<5 strong (k+a/a+k) 
galaxies in our models, we calculated the fraction of galaxies 
which have truncated their star formation shortly before the 
observed epoch: ie., t obs — t trU nc<1.5 Gyr (equivalent width 
Wo(H<5)>3 A) or t obs -t tIunc <l Gyr (W (R5)>5 A), and 
/ g as(£trunc)>0.1, where the calibration is based on Poggianti 
et al. (1999) and Balogh et al. (1999). These two definitions 
of the H5 strong galaxies correspond to the different defini- 



© 0000 RAS, MNRAS 000, 000-000 



14 Kodama & Bower 




Figure 10. The distribution of the star formation truncation time for the galaxies with My<— 19 at 2=0. The truncation epoch tlsLac 
is measured from the present-day rather than from the epoch at which the cluster is observed. The vertical axis gives the percentage 
of galaxies with this truncation time per Gyr. The bin sizes are set so that colour differences can be resolved within the photometric 
errors. Two cases are considered regarding the star formation at the time of galaxy accretion. In Panel (a), we assume that half of the 
remaining gas in the disk is lost out of the system and the other half is turned into stars in a burst at the galaxy accretion, and no 
residual star formation is present afterwards. In Panel (b), we assume that no gas is lost from the disk and only 1 per cent of disk gas 
is converted into stars by a burst. The remaining gas is consumed exponentially with an e-folding time of 1 Gyr. Thick dashed curves 
show the predicted galaxy accretion history of a 10 Mq cluster at z=0.33 for a CDM-like power in an open universe (f2o=0.2, <rg=0.8, 
r=0.1). The two curves in each plot are for threshold galaxy group mass of 5 ■ 10 12 Mq (right) and 5 ■ 10 13 Mq (left), respectively. 



Table 2. The fraction (%) of H<5 strong galaxies and [On] emis- 
sion line galaxies in the cluster cores. The MORPHS numbers and 
the CNOC numbers are taken from Poggianti et al. (1999) and 
Balogh et al. (1999), respectively. The 'corrected' CNOC numbers 
take into account the scatter due to measurment errors. 



obs/model 


H<5 (>3A) 


H<5 (>5A) 


[On] (>5A) 


MORPHS 


21 ± 2 


12 ± 1 


30 ± 5 


CNOC 


7.5 ± 1 


4.4 ± 0.7 


20 ± 5 


CNOC corrected 


< 2.5 


1.5 ± 0.8 




Model (a) 


6 


4 


ot 


Model (b) 


< 3 


< 1 


13 



Note - t If there is a time delay of 1 Gyr for the star burst 
epoch with respect to the accretion epoch, this fraction increases 
to 12 %.. 



tions adopted by these two papers. We also estimated the 
fraction of galaxies with [On] in emission, using Kennicutt's 
(1992b) calibration: 

SFRpVIoyr- 1 ] > 6.75 x 1(T 12 Wo (Oil), (4) 

Lb,® 

where Lb and Lb.q are the blue luminosity of a galaxy and 
the sun, respectively. The resulting fraction of H5 strong 
galaxies and emission galaxies are summarised in Table 2, 
where they can be compared with the fractions observed 
by the MORPHS (Poggianti et al. 1999) and the CNOC 
(Balogh et al. 1999) collaborations. For this comparison, we 



have transformed the numbers between the 3 and 5 A crite- 
rion using the same factor given by Balogh et al. (1999). A 
large difference in the fraction of H5 strong galaxies is evi- 
dent between the MORPHS and the CNOC numbers, with 
the MORPHS values being higher by factor 3. This differ- 
ence is still an open question, although the possibilities are 
extensively discussed in Balogh et al. (1999) including the 
intrinsic difference between the optical (MORPHS) and the 
X-ray (CNOC) selection of clusters. Balogh et al. also pre- 
sented the 'corrected' fraction of H5 strong galaxies, which 
has been corrected for the scattering effect due to the mea- 
surement errors. This effect reduces the fraction by about 
1/3. 

Our Model (b) gives numbers which are consistent with 
both the CNOC corrected fraction of H<5 strong galaxies and 
the fraction of [On] emission galaxies. This suggests that 
Model (b) is a good representation of the CNOC results even 
from spectroscopic point of view. Model (b) fails to repro- 
duce the fraction of H<5 strong galaxies found by MORPHS. 

In contrast, if we assume that the same correction for 
the scattering effect can be applied to the MORPHS num- 
bers, they are in good agreement with Model (a). However, 
Model (a) does not produce any [On] emission galaxies be- 
cause of the sharp truncation of star formation, hence no on- 
going star formation at the observed epoch. To reconcile this 
discrepancy, we allow for a time delay between the accretion 
of the galaxy into the cluster and the star-burst/truncation 
event. If the delay is 1 Gyr (the cluster dynamical time), 
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the fraction of [On] emission galaxies could be up to 12 per 
cent: still a factor of 2-3 lower than the observed fraction, 
but within the present uncertainties. 



5.4 The History of Star Formation in Cluster 
Galaxies 

Finally, we use our approach to compare the global star 
formation history in the cluster environment with that of the 
universe as a whole. The global star formation history of the 
universe has received extensive attention since Madau et al. 
(1996) combined data on the local Ha luminosity (Gallego et 
al. 1995), the UV luminosity density from the CFRS survey 
(Lilly et al. 1995) and the Lyman Break galaxies (Steidel 
et al. 1995). This paper is a first attempt to show a similar 
plot for the galaxies in the cluster environment. Although 
we do not have direct observations at z>0.5, the data still 
probe star formation histories out to z~l due to reasonably 
good sensitivity of the rest frame U—V colour. For each 
cluster, we sum up individual star formation histories over 
all the cluster members (except the brightest cluster galaxy) 
and normalised to a galaxy mass of 10 12 M Q . In Fig 11 we 
show the integrated star formation histories for the galaxies 
in the rich cluster cores within R30 for the four clusters 
(Coma, A2390, MS1358+62, and MS1621+26). Two panels 
correspond to the two extreme models (a: burst model, b: 
slow decay model) as in Fig. 10. 

The long-short dashed and the long dashed curves show 
the evolution of the field galaxy star formation as suggested 
by Madau et al (1998) and as fitted by our field galaxy 
model (Appendix B). We have assumed that the global star 
formation rate is constant beyond z=1.5. 

In the cluster cores, star formation at low redshifts de- 
clines much more rapidly than in the field, with galaxy for- 
mation skewed to higher redshifts. At 2 = 0.5, the typical 
star formation rate per galaxy mass is a factor 3-10 smaller 
in cluster cores, depending on the model chosen. This is in 
qualitative agreement with Balogh et al. (1998) and Couch 
ct al. (2000) who find a strong decline in the star formation 
rate from the field to the cluster environment at intermedi- 
ate redshifts. Although we prefer Model (b), which gives the 
most consistent picture of the assembling history of cluster 
galaxies, we should note that the exact locus of cluster star 
formation history on this plot depends on the treatment of 
the truncation, especially on the amount of gas that is lost 
from the system when it is accreted into the cluster (ie., 
eioss)- This uncertainty dominates the highest redshift bin 
since we have only a model dependent constraint (from the 
total number of galaxies now lying on the red sequence) on 
the accretion rate long before the cluster is observed. The 
important next step is to directly map the accretion history 
at higher redshifts by extending the method to more distant 
clusters. In addition, better understanding of the way in 
which star formation is truncated in the clusters and within 
individual galaxies eg., from Ha star formation surveys (eg., 
Balogh & Morris 2000; Couch et al. 2000) will better con- 
strain the modelling that is required to determine the cluster 
infall rate. 



6 CONCLUSIONS 

We have created the field corrected CM diagrams for the 

7 CNOC distant clusters (0.23£z<0.43) in order to study 
their galaxy populations in the same manner as in B084. 
We make a comparison with the redshift using data for the 
Coma cluster from Terlevich et al. (2000). The form for the 
redshift evolution of the blue galaxy fraction agrees well with 
the trend seen in B084. In addition, the blue galaxy fraction 
is generally a strong increasing function of distance from the 
cluster centre, and depends only weakly on magnitude. 

We have then applied these data-sets to build up a se- 
quence of snap-shots of cluster evolution. We follow the evo- 
lution of the whole population of cluster galaxies in these 
CM diagrams accounting for the flow of galaxies across the 
diagram as they are accreted by the cluster. 

• Firstly, we have shown that the blue galaxies are incor- 
porated into the present-day tight CM relation as they fade 
and become redder rapidly after the truncation of their star 
formation, giving the possible link to the faint (>M»+1) 
SO galaxies in present-day clusters. The brighter SO galax- 
ies, however, must have ceased star formation earlier than 
z ~ 0.4. In clusters, as in the universe as a whole, fainter 
galaxies have a more extended star formation history than 
their bright counterparts. Nevertheless, this effect does not 
lead to a sizable age variation along the CM relation at the 
present-day. 

• Secondly, we have used the distribution of galaxies in 
the colour space to infer the evolution of the rate at which 
galaxies have been accreted by the cluster. We interpret 
the observed distribution in colour by following the declin- 
ing star formation rate in the field and statistically assign- 
ing a star formation truncation epoch. If star formation is 
truncated abruptly, we find that too many recent accretion 
events are required relative to the past rate. In order to be 
consistent with a smoothly declining accretion rate, such as 
suggested by the extended P-S theory, we need to allow for 
some level of residual star formation after galaxy accretion. 
This suggests that ram- pressure and/or mergers/harassment 
cannot be so effective as to remove all the remaining gas 
from the system in a single orbital time-scale. The Butcher- 
Oemler effect is thus shown to result from a combination of 
the three effects, namely, declining accretion rate, declining 
field star formation rate, and the truncation of star forma- 
tion. 

• Finally, we have presented the global star formation 
history for the galaxies in rich cluster cores for the first time. 
It is shown that the star formation in cluster galaxies is much 
less than that in the field galaxies below z=l. However, the 
factor by which star formation is suppressed depends on the 
effectiveness of the cluster environment in suppressing star 
formation. 

It is important to ask how robust these conclusions are 
to the assumptions that must be made in order to model the 
galaxy properties. When propagating the observed proper- 
ties of galaxies forward in time, the past star formation rate 
has only a relatively weak effect (as shown by Fig. 6) and 
thus the results derived from this analysis are robust. De- 
termining the distribution of truncation times, and hence 
comparing the past and present accretion rates in clusters is 
significantly more model dependent. Firstly, the ratio of the 
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Figure 11. The star formation rate as a function of redshift for galaxies in rich cluster cores. The star formation rates have been derived 
from Model (a) (left panel) and Model (b) (right panel). The star formation rate (SFR) is presented per galaxy mass of 10 12 M Q . Four 
clusters are plotted individually by the curves with horizontal error bars showing the bin sizes; Coma (dot-dashed), A2390 (dotted), 
MS1358+62 (solid), and MS1621+26 (dashed). The long dashed line above indicates the star formation history for the CNOC field 
galaxies. The long-short dashed line above show the cosmic star formation history (field) taken from Madau ct al. (1998). 



present to past star formation rate must be estimated for 
the field galaxy population. This must be parameterised by 
a simple model and matched to the observed colour distribu- 
tion of field galaxies. The properties of the accreting galaxies 
can, however, be constrained by a wide range of observa- 
tions in addition to the observed field galaxy colours. The 
need to match the cosmic star formation history suggested 
by Madau et al. (1998) has motivated our choice of the gas 
infall model. Secondly, we must parameterise the effect of 
the cluster environment on the galaxies that are accreted. 
This is the more critical component of the model, and we 
illustrate this by considering two scenarios. If star formation 
is cut-off abruptly by a star burst and gas stripping upon 
accretion into the cluster, the model struggles to explain the 
high fraction of blue objects seen in the distant clusters. In 
contrast an extended period of star formation decay natu- 
rally reproduces a smooth decay in the accretion rate that 
is similar in all of the clusters. More complex scenarios arc 
easily possible: for example galaxies might either undergo a 
star burst or might decay slowly, depending on their orbit 
within the cluster. An important avenue for future work is to 
integrate spectral star formation indicators (such as the [On] 
or Ha line strength), or even radio continuum flux (Smail et 
al. 1999), and H<5 absorption line strength into this analy- 
sis so that the degeneracy between colour and past/present 
star formation rate can be broken. This will improve our 
understanding of the truncation process. 

So far we have not considered the effects of dust ob- 
scuration. This is a potentially important complication that 
is extremely difficult to model in a straightforward man- 
ner. Adding a uniform dust extinction that depends on the 



present star formation rate reddens the observed colours, 
causing us to underestimate star formation rates. However, 
both the field and some of the cluster galaxies will be affected 
so that the reddening cancels out to a first approximation. 
On the other hand, if localised dust extinction accompanies 
a star-burst on entry to the cluster environment as suggested 
by Poggianti et al. (1999) and Smail et al. (1999), we would 
recover an uneven accretion history (similar to the behaviour 
shown by Model (a)) since the 'intrinsic' (ie., dust corrected) 
colour distribution would include many more 'blue' galaxies. 

Clearly the potential of these techniques are best ex- 
ploited by observing more and more distant clusters. To re- 
solve the star formation history in clusters at higher red- 
shifts (z>l), we need a good homogeneous sample of clus- 
ters at higher redshifts (z>0. 5-1.0). The Rosat deep X-ray 
cluster survey (Rosati et al. 1998) and the coming new sur- 
veys by Chandra and Newton with high sensitivity cameras 
will provide an ideal sample of high redshift clusters for this 
purpose. 

The analysis presented in this paper is limited by the 
accuracy of the the distant cluster photometry, as well as 
by the number of clusters in our sample. A larger sample 
of clusters at each redshift would allow projection effects to 
be averaged over and the variation in evolution history to 
be investigated. Wider spatial coverage is also important to 
more accurately subtract field population and also to inves- 
tigate the environmental effect beyond the cluster core out 
to the virial radius. This is essential if we are to improve 
our knowledge of the star formation truncation mechanism 
itself. At present, our approach is also compromised because 
we do not allow for the dynamical evolution of the cluster 
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populations. Because of this, we cannot compare the proper- 
ties of galaxies at different radii within the clusters. Recent 
rapid progress in the high resolution N-body simulation of 
clusters (eg., Moore et al. 1999b; Balogh et al. 2000; Diaferio 
et al. 2000) will allow this aspect to be taken into account. 
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APPENDIX A: STATISTICAL FIELD GALAXY 
SUBTRACTION 

This appendix describes our method for subtracting the fore- 
ground/background populations from the cluster fields. The 
sample of field galaxies is taken from the edge of the cluster 
fields of three CNOC clusters with wide imaging fields; ie., 



11 < ARA < 21 
8 < ADEC < 11.5 
and -12 < ADEC < -i 



for A2390, 

for MS1358+62, 

for MS1621+26, 



where ARA and ADEC are the distances from the cluster 
centre in the unit of arcmin, and these correspond to ^3 Mpc 
which is comparable to or slightly smaller than the virial ra- 
dius (Carlberg et al. 1996). We omit the MS1512+36 cluster 
because of probable contamination by a foreground system. 
We can estimate the contamination of the 'field' sample ei- 
ther by extrapolating the surface density profile of clusters 
or directly from the CNOC spectroscopic data. Both meth- 
ods suggest that only 10-15 per cent of the galaxies are 
expected to be associated with the outer regions of the clus- 
ters. After rejecting a few spectroscopically confirmed clus- 
ter members from this field area, the total number of the 
field galaxies is 871 down to r=23.5. Divided by the to- 
tal area of 137.5 arcmin 2 , it gives the field number density 
of 6.34±0.21 galaxies per arcmin . This can be compared 
to the deep field galaxy counts in Gunn's system (Brain- 
erd et al. 1995) which has 90.1 arcmin 2 field of view and 
7.01±0.28 galaxies per arcmin 2 with the same magnitude 
cut. Furthermore, both the magnitude and colour distribu- 
tion of galaxies are similar between the two fields. Therefore 
our field sample can be regarded as a good representation of 
the field population. Even if we used Brainerd et al.'s (1995) 
field population, the results presented in this paper would 
hardly change. The CNOC field galaxies defined in this way 
are plotted on the CM diagram in the bottom right panel of 
Fig. 2. 

We build the two dimensional distribution histogram of 
these field galaxies with bins of width 0.3 in colour and 0.5 
in magnitude. A similar histogram is also created for each 
cluster field within 7?3o. This obviously includes the field 
galaxies that are to be subtracted. The field histogram is 
normalised for each cluster to match the area within Rzo- 
The number of spectroscopically confirmed non-members is 
then subtracted from these two histograms in each bin. The 
resulting number can be negative. We now count the number 
of galaxies in each bin for the field (A^ eld ) and the clus- 
ter plus field (A^iustor+Hcid) The formcr can be larger than 

the latter for low density bins due to low number statistics. 
If this happens, we re-distribute the excess number of field 



galaxies to the neighbouring bins (8, 5, or 3 bins depending 
on the location of the bin under concern) with equal weight 
until 



Arficld a rclustcr+ficld 

ij ^ ij 



(Al) 



is satisfied in all bins. We then define a probability, P, for 
each galaxy in the cluster field that it is a field galaxy: 



V 



clustcr + ficld ' 



(A2) 



where the numbers are taken from the bin that the galaxy 
belongs to. We now apply Monte-Carlo simulations to sta- 
tistically determine the field galaxies. We randomly pick out 
a membership-unknown galaxy from the entire CM diagram 
and refer to its probability P to determine whether it should 
be regarded as a field galaxy. We repeat this process for all 
galaxies in the cluster field and rerun the Monte-Carlo sim- 
ulation to average over the 100 realisations for each cluster. 



APPENDIX B: THE INFALL MODEL 

We use the gas infall model prescription for star formation 
in galaxies under the instantaneous recycling approximation 
(Tinsley 1980). We do not take into account the chemical 
evolution in the galaxy system in this paper, and the metal- 
licity is fixed at solar value. This treatment is supported by 
BKT98, where they found the colour changing rate is only 
a weak function of metallicity. 

The galactic gas is supplied from a surrounding gas 
reservoir (halo) at the rate: 



fcn(t) = -^-exp (-±-\ 

'in V (in/ 



(Bl) 



The total mass of the gas in the reservoir is normalised to 
unity. The zero-point of the time t is set to the look back 
time of 15 Gyr which corresponds to the redshift of 5.4 in 
our cosmology. 

A Schmidt law (Schmidt 1959) with a power unity is 
adopted, where star formation rate is proportional to the 
gas density, hence mass, as 

m = -M gas (t), (B2) 

T * 

where M gas indicates the mass of the gas that has infallen 
and has not yet been turned into stars. The initial gas mass, 
Mgas(0), is set to zero. Note that this formulation gives an 
exponentially decaying star formation rate with an effective 
time scale r° ff = t,/q in the case of the simple models, 
where a is the so-called locked-up mass fraction defined by 
Tinsley (1980). The Salpeter mass function (Salpeter 1955) 
that we adopted in this paper with lower mass cut-off 0.1 Mq 
and upper mass cut-off 60 M , gives a — 0.72. 

In order that the gas fraction is stable, we adopt 
eff ... _ 



(B3) 



This assumption leads to the simple analytic solutions for 
the star formation in a galaxy: 

t' 

a V t / 



tp(t) = i-^exp 



(B4) 



which has a peak at t = r. 

The gas mass and stellar mass are given by: 
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M gas (i) = ^exp 



(B5) 



and 

M, 



.(t) = l-(l + |)exp(-|), (B6) 

respectively. We also define the gas mass fraction as: 

M gas (t) 



/gas (t) 



M gas (i) + M„(t)' 



(B7) 



Equations from (Bl) to (B6), except (B3), are scaled 
for each galaxy to match its absolute magnitude. 



APPENDIX C: 
BURST 



TRUNCATION AND STAR 



When a galaxy is accreted from the field environment by 
the cluster (£ a cc), all the gas in the halo is assumed to be 
stripped off due to tidal stripping and the ram-pressure of 
the intra-cluster medium. This prevents further gas being 
supplied from the halo gas reservoir. A substantial fraction 
of the gas in a galaxy disk (ie., the material which has been 
accreted from the halo) may also be stripped away by the 
ram-pressure as numerical simulations suggest (Abadi et al. 
1999; Quilis et al. 2000). The mass of gas that is lost from 
the system is described using the remaining gas mass and 
the gas loss efficiency parameter ei oss as 



Mo 



M gas (i acc J X <3lo 



(Cl) 



We allow a star burst to accompany the gas stripping 
process. Gas mass that is transformed into stars in the burst 
is described by a burst efficiency parameter e bu rst and is 
given by 



M bu 



(C2) 



The burst strength is represented by the mass fraction of 
the burst population to the total stellar mass: 

Afburst 



/"burst 



M„(tacc) +M burst ' 



(C3) 



After these processes, star formation in the galaxy is 
rapidly truncated at time itrunc. We assume that the star for- 
mation truncation time itrunc is approximately equal to the 
galaxy accretion time i acc - After itrunc, the galaxy quickly 
consumes the remaining gas in the disk following a Schmidt 
law: 

i>{t) = 1 ~ ei ° SS ~ SbUISt M gas (itrunc) CXp f J-J^ ) (04) 
£^T~trunc \ 



Ttr 



where rtrunc is a time scale for the truncation of star for- 
mation and is fixed to 1 Gyr. In principle the model could 
be adjusted to use a range of decay timescales. We adopt 
1 Gyr since it is midway between the average timescale for 
field galaxies and abrupt truncation. 

We consider two specific models in this paper that are 
representative of the extreme possibilities: 



Model 

Gloss 

(a) 0.5 

(b) 0.0 



^burst 

0.5 
0.01 



Model (a) is motivated by Quilis et al. (2000) and reflects 
what would happen if the diffuse hydrogen were removed 




Figure Dl. Distribution of g—r colour for the spectroscopically 
confirmed field galaxies which arc located between 0.3<2<0.35. 
These galaxies are taken from all of our 7 CNOC cluster fields 
and the total number amounts to 94. The dashed line corresponds 
to the Gaussian distribution of r for the field galaxies in equa- 
tion (Dl). 



from the galaxy by ram-pressure stripping while the remain- 
ing molecular gas (typically 50 % of the total gas mass) is 
consumed in a pressure induced star burst. Model (b) mim- 
ics the decline in star formation rate expected if the action 
of the cluster is simply to remove the outer gas reservoir. 



APPENDIX D: ASSIGNING THE STAR 
FORMATION TIME SCALE AND THE 
TRUNCATION EPOCH 

This section is devoted to explaining how we assign the star 
formation/gas infall time scale (r) and the truncation epoch 
(itrunc) for each cluster galaxy. 

Since we consider that the cluster galaxies is accreted 
from the field, the distribution of r for the cluster galaxies 
is taken to be the same as for the field population. This is 
directly estimated from the colour distribution of the real 
field galaxies through our treatment of star formation de- 
scribed in Appendix B. Here we use 94 spectroscopically 
confirmed field galaxies between 0.3<z<0.35 from the entire 
observed fields of the seven CNOC clusters. The estimated 
distribution of r for the field galaxies can be expressed by a 
Gaussian: 



$(t) oc cxp < — 



H( 



log r - log r * 



(Dl) 



where log to=0.3, and tr=0.2. This distribution not only 
gives a reasonable match to the colour distribution of the 
CNOC field galaxies (Fig. Dl), but also reproduces the so- 
called Madau curve (Madau et al. 1998); ie. the global star 
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formation rate for field galaxies as a function of redshift at 
the same time (see the long dashed curve against the long- 
short dashed curve in Fig. 11). 

When the galaxies fall into clusters, star formation 
ceases with either gas loss and/or a star burst as described in 
Appendix C. Now, in order to reproduce a certain colour of 
a given cluster galaxy, various combinations of r and itrunc 
are possible for a given set of ei oss and e b urst, although r 
should be larger than a certain minimum value r m i n which 
corresponds to ttrunc = 0. 

For each possible combination of (T,ttrunc), the colour 
changing rate R at the cluster redshift is calculated: 

Tt(rt \- d( - 9 ~ r) 

il[T, ItruncJ - 



dt 



(D2) 

- z clustcr 

Note the exception that the colour used is U—V for Coma. 
The probability of assigning t, q(r), is inversely proportional 
to this rate since a galaxy with higher colour changing rate 
is rarer to find: 

q( r ) — J C /^( T i 'trunc) T > T~min (D3) 
Hy ' 1 T < 

where the normalisation c is determined by 

q(r)d\ogT = 1. (D4) 



The distribution of r is thus skewed towards larger val- 
ues due to the minimum cut r m i n and the q(r) term. To get 
the same distribution of r in Equation (Dl) for the clus- 
ter galaxies, the probability of assigning a given r should 
be corrected for these effects. The correction function Q(t) 
is the sum of the individual q(r) term over all the cluster 
members: 

Q(r) = E 9i (r). (D5) 

Finally, we assign r randomly for each cluster galaxy 
using the following distribution function $'(r): 

$'(r) = $(r)/Q(r) X q(r). (D6) 

The truncation time, ttrunc, is then determined by requiring 
the galaxy to have its observed colour at the cluster look- 
back time. 
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